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We study the free energy distribution function of weakly disordered Ising ferromagnet in terms 
of the 73-dimensional random temperature Ginzburg-Landau Hamiltonian. It is shown that besides 
the usual Gaussian "body" this distribution function exhibits non-Gaussian tails both in the param- 
agnetic and in the ferromagnetic phases. Explicit asymptotic expressions for these tails are derived. 
It is demonstrated that the tails are strongly asymmetric: the left tail (for large negative values of 
the free energy) is much more slow than the right one (for large positive values of the free energy). 
It is argued that in the critical point the free energy of the random Ising ferromagnet in dimensions 
D < 4 is described by a non-trivial universal distribution function being non self-averaging. 
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I. INTRODUCTION 



The ferromagnetic systems with weak disorder (e.g. with small concentration of non-magnetic impurities) are under 
intensive study during last four decades. As the disorder of this type can not destroy the ferromagnetic ground state, 
most of the efforts (both theoretical, experimental and numerical) were concentrated on the critical properties of 
such systems (see e.g. [l|, 0]) where even weak disorder was proved to be "relevant variable". It turns out that 
the presence of quenched disorder can essentially modify the critical behavior of the system such that new universal 
critical exponents may set in. Another interesting phenomena which appear due to quenched randomness in the Ising 
ferromagnet are the so called Griffiths singularities in the thermodynamical functions 0. Apart from these two big 
issues the effects produced by weak quenched disorder on thermodynamical properties of a ferromagnet (away from 
the critical point) were usually considered as more or less irrelevant corrections to the pure case. 

The aim of the present study is to demonstrate that due to the presence of weak disorder the statistical properties 
of (random) free energy of the Ising ferromagnet are getting rather non-trivial. Of course, in the thermodynamic 
limit due to self-averaging the free energy distribution function turns into the (5-function. However, if the volume of 
the system is large but finite (which, in particular, is the case in numerical simulations) the situation is much more 
complicated. Away from the critical point at scales much bigger than the correlation length Rc the system could be 
considered as a set of more or less independent regions with the size Rc- In this case one would naively expect that 
the free energy distribution function must be Gaussian. In fact this is not the case. It can be shown that besides the 
central Gaussian part this distribution has asymmetric and essentially non-Gaussian tails. Approaching the critical 
point one finds that the range of validity of the Gaussian "body" shrinks while the "tails" are getting of the same order 
as the "body". Moreover, if we naively imitate the critical point by assuming that the correlation length becomes of 
the order of the system size we formally find that the free energy distribution function turns into a universal curve. 

In our present study the Ising ferromagnet with quenched disorder will be considered in terms of the random 
temperature Z?-dimensional Ginzburg-Landau Hamiltonian: 

= j d^x[i(V</.(x))Vi(r-^(x))02(x) + i5o</''(x)] (1) 

where independent random quenched parameters ^(x) are described by the Gaussian distribution function 

1 / 1 



pu^)] = n 



(2) 



where Uo is the parameter which describes the strength of the disorder, = 2mo. In what follows it will be supposed 
that the disorder is weak: Uo ^ and Uo go- 

For a given realization of the disorder the partition function of the considered system is 

Z[^] ^ Jt^<P exp(~i?[</),e]) - exp(-^^[C]) (3) 

where J Dcf) denotes the integration over all configurations of the fields (/'(x) and is the free energy of the system. 

The distribution function of the random quantity F [^] can be analyzed by studying the moments of the partition 
function. Taking the integer n-th power of the expression in eq.([3|) and averaging over the disorder parameters ^(x) 
(with the Gaussian distribution, eq.([2])) we get the replica partition function 

=Z{n) = j Vcpi... j V4>,,eiq>{~H^[ct>\) (4) 

where the replica Hamiltonian 

i/„[0] = /'d^x[iX^(V0,(x))Virf^</.2(x) + i 5a6</>^(x)0?(x)] (5) 

a=l a=l a, 6=1 

depends on n interacting fields 4> = {^ij ■■■4>ri\, and the coupling {n x n) matrix 



gab = goSab - Uo 



(6) 
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Taking into account the definition of the free energy in eq.([3]), the rephca partition function, eq.(l4]), can also be 
represented as follows, 



Z{n) = exp(-nF[C]) (7) 

or 

/ + 00 
dF P{F) eyip{-nF) (8) 
-oo 

where P{F) is a free energy distribution function. 

Usually the replica partition function, eq.Q, is studied for deriving the average value of the free energy. The 
heuristic (not well justified) procedure of the replica calculations requires to perform the analytic continuation of the 
function Z{n) from integer to arbitrary values of the replica parameter n. Then, taking the limit n — >■ in both sides 
of eq.® (and taking into account that J dFP{F) — 1) we formally get 

liin^M_i.^ r dFP{F)F = F (9) 

In fact if we are interested not just in the average free energy of the system but in the properties of its distribution 
function the situation even with weak disorder turns out to be rather nontrivial. Considering the system away from 
the critical point one would naively expect that at least at scales greater than the correlation length (where the system 
is expected to split into a set of more or less independent "cells" of the size of the correlation length) "everything 
must be Gaussian distributed" . One can easily prove that this is not true. Indeed, let us consider the coupling matrix 
gab, eq®. It has (n — 1) eigenvalues Ai = go — Uq and one eigenvalue X2 = go — nuo- We see that for n > go/uo the 
eigenvalue A2 is getting negative. Thus, according to eqs.(|4])-([5|), we find that all moments of the partition function 
Z{n) with n > go/uo must be divergent! On the other hand, according to the relation this indicate that the free 
energy distribution function of the considered system can not be Gaussian. Moreover, since it is evident that the 
divergency in the integration in the rhs of eq.® takes place in the limit F — ^ — cxi, we can easily guess that the left 
asymptotic of the distribution function P{F) has the following simple form: 

P(F^ -cx)) - exp(-— IfI) (10) 



It should be stressed that this phenomenon must be quite general: it take place for any dimension of the system 
independently of the value of the temperature parameter t (except for the critical point) both in paramagnetic and 
in the ferromagnetic phases. 



II. TOY MODEL 



To analyze the properties of the free energy distribution function of a random ferromagnet one can start with very 
simple "toy" model. Let us consider the system described by the Hamiltonian containing only one degree of freedom: 

H{cl>,0 = l{r^Oc^' + \gocj^' (11) 
The random quenched parameter ^ is described by the Gaussian distribution function 




and its typical value = 2uo is supposed to be small, Uo ^ t'^ and Uq <C go- Besides, to simplify the calculations it 
will be assumed that the coupling parameter go is also small, go ^ ■ 



A. "Paramagnetic" region (r > 0) 



To simplify notations let us redefine: ^ — > t^, Uo 
partition function of the system is 



T^w, go T^g and (j) — > (^/^/r. It these new notations the 



Zparai^^') 



+00 



exp 



exp [-FparaiCl] 



(13) 
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where FparaiO (random) free energy of the system. For the rephca partition function we get: 



+00 



'Anu 



(14) 



On the other hand, 



/+cx> 
dFPparaiP) exp[-7li^] 
-00 



(15) 



where Ppara (F) is a free energy distribution function. The above relation is the bilateral Laplace transform. Performing 
analytic continuation of the function Zpara{n) from integer n to arbitrary complex values, by inverse Laplace transform 
we obtain: 



Ppara{F) = / — Zpara{z) exp[z F] 

J -ioo ^'^^ 



(16) 



Under conditions u ^ g ^ 1 one can easily perform (approximate) calculations of the partition function, eq.(|13p. For 
^ < 1 the main contribution in the integration over comes from the vicinity of the point = 0. On the other hand, 
at ^ > 1 the integral is dominated by the vicinity of the "ferromagnetic" points (j) = ±-^/(^^^^TJ73. Thus, 



Zp^^aiO expf-/o(0 + ^^-r^e(e-i; 



where 



/o(0 



ilog|e-l|, for V5 
ilog(5), for Vff 



(17) 



(18) 



Correspondingly, for the rephca partition function, eq. (|14p . we find 



Zpara (j^^ 



+00 



' Attu 



exp 



1 

Au 



e-nFparaiO 



where 



FparaiO ^ MO ~ ^^^0{^ - 1) 



(19) 



(20) 



One can easily note here that since Fpara{£, +00) ^ ^J^^^ ^^'^ partition function Zpara{n) in eg. (1191) is divergent 
for all values n > g/u. 

Finally, substituting eqs. (|T9|) - (l20|) into eq. (|16l) (where the integration contour can be chosen to be just the imaginary 
axes, z = iio) we obtain 



Ppara {F^ 



-00 2^ 
+°° d^ 



-|-oo 



/ Attu 
1 



exp 



+iuj{F~Fpara{0) 



/ Attu 

>(+) 



: exp 



Au 



S[F~FparaiO) 



(21) 



Right part of this distribution function Ppara{F) valid for — ilog(l/.g) ^ F < +00 is defined by the region 
-00 < (^ - 1) < where FparaiO ^ \ log(l - C)- Using eq.(in|) we get 



p(+) (p) 

para\ / 



exp 



In the vicinity of the point F^ = u the above distribution has the Gaussian form: 

1 



p(G) (p) 

para\ J 



exp 



^{F-uY 



(22) 



(23) 
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However, the above form of the distribution is vahd only for |F| <C 1. For large (positive) values of the free energy 
the distribution becomes essentially non-Gaussian: 



Ppara{F ^ +0o) 



exp 



(24) 



which goes to zero much faster than the Gaussian curve. 

Left asymptotic of the distribution function, Ppara{F) valid for — oo F <S^ —j log(l/5) is defined by the region 
V5 < (C - 1) < where Fpara{0 ~ -:^(e - 1)^- Using eq.(l2T]) we get: 



(F) 

para \ J 



9 



exp 



1 



47rM|F| 

Thus we find that the left tail of the free energy distribution function Ppara{F) is also non-Gaussian: 

Ppara{F^-(^) ^ exp(~^|F|) 



(25) 



(26) 



and it goes to zero much slower that the Gaussian curve. 

In conclusion we find that the free energy distribution function of the "toy" model, ea.pTj). at r > is essentially 
non-Gaussian and strongly asymmetric: its right tail, eq. (p4|) . is much steeper than the left one, ea. (l26l) . The example 
of the entire distribution function Ppara{F) for g = 0.1 and u — 0.01 is shown in Figure 1. 




FIG. 1. Free energy distribution function of the "paramagnetic" toy model, ea. (fT3|l . with g = 0.1 and u = 0.01 



B. "Ferromagnetic" region (r < 0) 

For the negative values of the temperature parameter r in the toy model, eg. ([TT]) . we redefine: ^ 



go ^ T g and 



t|. It these new notations the partition function of the system becomes 



exp 



exp[-F/erro(C)] 



(27) 
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where random parameters ^ are described by the Gaussian distribution, eq.([T2]). The calculations for the "ferromag- 
netic" case are quite similar to the case r > 0, considered in the previous subsection. In particular, for the partition 
function we get: 



+00 



where 



Zferro{n) ~ / CXp - - nF/erro (C) (28) 



Fj^rroiO ^ h{O-^^-^0{l-O (29) 



Here /o(^) is defined in eg. ([T8| . As in the "paramagnetic" case, one can easily note that since ^/erro(f ^ — oo) 
— ^^/4(7, the partition function Zferro{n), eq.(|28p. with n > g/u is divergent. 
Similarly to ea.(PT|). for the free energy distribution function we get: 



+00 



-—^■'+tUj{F-Fferro(.0) 



+ 00 



/ Attu 



exp 



6{F-F, 



ferro 



(30) 



The right tail of this function, Pftrroi^)^ vahd for —j log(l/.g) <^ F < +00 is defined by the region ^/g <^ (C^l) < +00 
where -F/erro(C) — 5 log(^ ~ !)• Using eg. ([30)) we find 



Pferro i-^) 



exp 



2F- -3-fexp(2i^) +1 



(31) 



Thus the decay of the right tail of this function coincides with the one of the "paramagnetic" function, ea ([M|) . 



Pferro{F ^ +00) ~ CXp -— exp(4F) 



(32) 



The left tail of the free energy distribution function, ^ye„o(-^)' valid for —00 < F — ilog(l/5) is defined by the 
region —00 < (^ — 1) ^ ^Vd where Fferro{£.) — ~ j"(l ^ C)^- Using eg. ([50]) we find 



P. 



ferro 



exp 



1 

iu 



4ttu\F\ 

In the vicinity of the " ferromagnetic" average free energy F^ ~ — 1/4(7 this distribution is Gaussian, 

"2 1.21 



P 



(G) 



ferro 



(F) 



exp 



F 



(33) 



(34) 



However at large negative values of the free energy (similarly to the " paramagnetic" case) the distribution becomes 
essentially non-Gaussian: 



Pferro{F^-^) - CXp(-^|F|) 



(35) 



The example of the entire distribution function Pferro{F) for g = 0.1 and u ~ 0.01 is shown in Figure 2. Again, like 
in the "paramagnetic" case this function is strongly asymmetric: its left tail, eg. (|35|) . is much slower that the right 
one, eg. (|32)) . 



III. MACROSCOPIC SYSTEM 



Let us consider _D-dimensional random ferromagnet described by the Ginsburg-Landau Hamiltonian, eg.(IT]). In 
the mean-field approximation away from the critical point (in the paramagnetic phase, r > 0) this system can be 
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FIG. 2. Free energy distribution function of the "ferromagnetic" toy model, ea. (l27|l . with g = 0.1 and u = 0.01 



considered as a set of ~ V/R^{t) independent "sells", where V is the volume of the system and Rc{t) ~ 1/-\/t 
is the correlation length. Performing spatial rescaling x — >■ RcX = T~^/^a;, according to the standard ideas of the 
scaling theory we have to rescale (j> — >■ t^^~^^/^(/<, go — T~''^~^^/^g and Uo — )■ t~^^~''^/^u. Then in terms of the new 
(rescaled) parameters the partition function of the considered system can be estimated as follows 



N 



m - n 



d4) 

J2^ 



exp 



N 



(36) 



k=l 



where the function depends on N random parameters ^ — {^i, ...^tv} and Zpara{C) is the partition function of 

the toy model, ea.(|13p. The parameters {f/t} are described by the independent Gaussian distributions 



N 



p[^\ = n 



k=l 



1 

Au' 



(37) 



Correspondingly, for the replica partition function we obtain 



Z(n) 



4-00 



-1 AT ^ 



where Zpara{n) is the replica partition function of the toy model, cq. ()14p . On the other hand. 



/ + OC 
dFV{F) exp{-nF) 
-OO 



(38) 



(39) 



where ViF) is the total free energy distribution function of the considered macroscopic system. Performing inverse 
Laplace transform we get 



V{F) ^ / IZparaiz)) exp{zF) 



(40) 



8 



Taking z = iw, introducing the free energy density, / = F/N, and substituting here ea. (|19p . for the corresponding 
probabihty distribution function we find the following expression 



f+°° du \ /■+°° d^ / 1 .9 r . 



N 



(41) 



where the function FparaiO is defined in eqs.ipUl and 

The main Gaussian part of the above distribution valid for the values |/| ^ 1 is defined by the region |^| ^ 1 where 
FparaiO - \ log(l - C) — "H" ^his case 



+^ dw 



.^=exp + ■7i^^ + ^^f 

/Anu V 4u 2 



N 



1 



icxp 



Correspondingly, for the distribution function of the total free energy F = JN — fV/R^ we get: 



RP 



If 

Vttu ^ V Vu 



exp 



i-^F^) 
\ Vu J 



(42) 



(43) 



One can easily note that in the thermodynamic limit — > oo this distribution function turns into the (5-function 
which assumes that the considered random system is self-averaging. However, at large but finite volume of the system 
one can easily prove that both the left and the right tails of the function 'P{F) are essentially non-Gaussian. 

Right tail V^^\f) of the distribution function, ea.(|4T|). valid for the values / 1 is defined by the region |^| ^ 1 
(^ < 0) where Fpara{£,) — \ log |^|. In this case 



N 



75(+)(/)^n 



fe=i 



exp 



+ 00 



/Attu 



exp 



k=l k=l 



N 
4m 



exp (4/) 



exp 



V 



AuRi 



■exp (4/) 



Correspondingly, for the right tail of the total free energy distribution function we get 



p(+)(F) ~ exp 



V 



{4") 



(44) 



(45) 



which is valid for F ^ V/R^. 

Left tail V^^\f) of the distribution function, ea.(PT|). valid for / < and |/| 3> 1 is defined by the region ^ 3> 1 

par 



where FparaiO - ^T^(^- ^^^^ case 



+'°° dz_ 

-i^ '2,711 



+00 



d^ 



+'°° dz_ 



'4:7TU 

-N/2 



exp 



4wV g 



N 



exp{Nf2 



l--zj exp(-iV|/|z) 



(46) 



Shifting here the integration contour to the right by the value g/u: z — g / u + iuj and neglecting pre-exponential 
factor, we obtain 



V<^-\f) ^ exp(-7V%|^ 
or, for the total free energy (the result already discussed in the Introduction, eg. dTUl) ). 

V'^^Hf) ^ exp(-Jli^l) 

which is valid for |F| > . 



(47) 



(48) 
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Thus, according to eas. (|43|) . (|45l) and (l48l) for the entire free energy distribution function of the random Ising 
ferromagnet (in the paramagnetic phase away from the critical point) we find the following result: 



'exp(-flFl); 



for F < and |F| > 



r{F) 



exp I 



exp 



Vu- 



W^xp 



for |F| « ^ 



for F > 



(49) 



where (in the framework of the mean-field scaling theory) Rc — l/Vr is the correlation length, V is the volume of the 
system, u = t^-^^^'^/^mq is the rescaled strength of the disorder, eq.([2]), and g = t'^^^'^'/^^q is the rescaled coupling 
constant of the original Ginsburg-Landau Hamiltonian, eq.([T|). 

One can easily repeat similar calculations for the free energy distribution function in the ferromagnetic phase to 
obtain the result similar to eq. p9)) . The only difference is that in this case the Gaussian central part of the distribution 
will be concentrated around the ferromagnetic ground state energy Fq — —^{y/R^) = ^t^^- 



IV. CONCLUSIONS 



One can easily note at least three apparent consequences of the obtained result, eg. (|49| . for the free energy distri- 
bution function of the random Ising ferromagnet. 

First. In the thermodynamic limit, — > oo, this distribution turns into the 5- function. This is not surprising and 
it just indicates that the considered system is self-averaging. 

Second. At large but finite volume of the system in addition to the trivial Gaussian "body" the distribution function 
exhibits essentially non-Gaussian tails. Moreover, these tails are strongly asymmetric: the left tail, at F — —oo, is 
much slower that the right one, at F -> +oo. It is this type of asymmetry of the free energv distribution function (the 
left tail is slow while the right tail is fast) which is observed in others random systems It is interesting to note 
that the result, ea.(|49)). has to be very general taking place at all temperatures (except for the critical region) and 
in all dimensions. On the other hand the behavior of each tail (both left and right) looks rather universal. Indeed, 
at large negative values of the free energy the logarithm of the distribution function is expected to be linear in |F|, 
while at large positive energies, it is the double logarithm of the distribution function which is expected to be linear 
in F. Presumably these two types of asymptotics could be relatively easy recovered by computer simulations. 

Third. While approaching the critical point we observe that the range of validity of the Gaussian "body" of the 
distribution function shrinks while the "tails" are getting of the same order as the "body". Of course, simple mean- 
field calculations considered in this paper are not valid in the critical point. Nevertheless, if we naively imitate the 
critical point by assuming that the correlation length becomes of the order of the system size (i.e. R^ ~ V) while the 
coupling parameters g and u (in dimension D < 4) become equal to their universal (depending only of the dimension 
D) fixed point values and u*, [l[ we find that (in the thermodynamic limit V — >■ oo) the free energy distribution 
function turns into the universal curve such that 



' exp(-f^|F|); forF < Oand |F| > 1 

(-^F2); for|F|«l 



exp 



exp 



(50) 



■4^exp(4F)^ 



for F » 1 



All that evokes the old idea that free energy (as well as some others thermodynamic quantities) of the random 
ferromagnets could be non self-averaging in the critical point @. Of course, the systematic derivation of the free 
energy distribution function in the critical point requires a special consideration. 
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